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Abstract 

We investigate quantum chromodynamics in 2+1 dimensions (QCD 3 ) 
using the Hamiltonian lattice held theory approach. The long wavelength 
structure of the ground state, which is closely related to the confinement phe¬ 
nomenon, is analyzed and its vacuum wave function is evaluated by means 
of the recently developed truncated eigenvalue equation method. The third 
order estimations show nice scaling for the physical quantities. 
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QCD has been accepted to be the most successful gauge theory of strongly 
interacting particles. In QCD, hadronic matter is composed of quarks, and 
interactions between them are mediated by eight massless gluons generated 
by the SU(3) gauge group. Asymptotic freedom of QCD at short distances 
makes the perturbative calculations of high energy processes possible. At 
long distances, however, there are a lot of low energy fundamental properties 
like confinement of quarks and gluons, vacuum structure, chiral-symmetry 
breaking, glueball masses, hadronic spectrum and weak interaction processes, 
and behaviors of hadronic matters at high temperature or high density, which 
can not be studied perturbatively. 

Over the past two decades, lattice gauge theory (LGT) has developed into 
a promising first principle nonperturbative approach to these phenomena. 
QCD in the pure gauge sector possesses a nontrivial vacuum structure and 
bound states called glueballs. The Hamiltonian LGT provides a convenient 
tool for the estimations of the wave functions for the ground state and excited 
states. Furthermore, physical observables such as glueball masses correspond 
to the eigenvalues. In this context, both numerical [1, 2, 3] and analytical 
[4, 5, 6, 7, 8, 9] efforts have been made. In [3], Arisue argues that the wave 
function of the ground state of a D-dimensional non-abelian theory for the 
long wavelength configurations A is 

TM) = exp[— f d D ~ 1 x trF 2 — ^ f d D ~ 1 x tr(VF) 2 ] } (1) 

e 2 J e b J 

with e being the renormalized coupling, F the held strength tensor and T> the 
covariant derivative. The idea was confirmed by his Monte Carlo simulation 
of a (2+l)-dimensional SU(2) lattice model. As well explained in [3], the 
correlation length of the continuum gauge held strength tensor in the vacuum, 
given by the square root of the ratio of the coefficients in (1), has the order 
of 1/e 2 , the same order of the conhnement scale of the theory. 

The motivations for the investigations of the (2+l)-dimensional lattice 
models are as follows. 

(1) they have potential applications to high temperature superconductivity; 

(2) they have many similarities to QCD in 3+1 dimensions (QCD 4 ) like 
asymptotic freedom, quark conhnement, spontaneous chiral-symmetry brea¬ 
king, and meson and glueball spectrum. It would be more economical to 
test various techniques on such a nontrivial theory which has as much as 
the same properties as QCD 4 but simpler (due to the advantage of lower 
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dimensionality and superrenormalizability). 

Recently, we proposed an analytical method [6, 7] for understanding the 
long wavelength behavior of Hamiltonian LGT, which is similar to Greensite’s 
method [4] of truncated eigenvalue equation but with a different truncation 
scheme. The vacuum wave function of (2+l)-dimensional SU(2) gauge theory 
were evaluated up to the third and fourth orders, and nice scaling behavior 
and agreement with the Monte Carlo data [3] were observed. 

It is very desirable to extend the method to a more realistic gauge group: 
SU(3). As will be seen, because of the nature of the group, the classification 
of the graphs is more complicated. To our knowledge, there has not been 
a published work on detailed investigation of the long wavelength vacuum 
structure of pure QCD 3 at zero temperature. There exist only preliminary 
analytic calculations in QCD 3 with fermions [10], in addition to some analytic 
analysis in the continuum [11] and Monte Carlo data [12, 13] for quenched 
QCD 3 at finite temperature. Recently, there has been an attempt [14] to 
include dynamical fermions in the numerical simulation of QCD 3 . 

The purpose of this paper is to describe and further explore our method 
for studying the vacuum structure of QCD 3 in the pure gauge sector. The 
study of the ground state properties is a first step towards the understanding 
of the structure of the glueballs and hadrons. The starting point is the dis¬ 
cretization of the Yang-Mills theory in the Hamiltonian formulation (discrete 
in space, continuous in time, and temporal gauge A(x, k t ) = 0) 

H = f £ E ? E ? - £ Tr(U r + U' p - 2), (2) 

la i ag z ^ 

where g is the bare coupling, a is the lattice spacing, Ef = E a (x,k ) is 
the color-electric held on the link / at site x and positive direction k, and 
the second term is the color-magnetic energy with U p being the product of 
link variables Ui = exp[igaA(x } k)\ around an elementary plaquette. In 2+1 
dimensions, the bare coupling and the lattice spacing have a simple relation 
g 2 = e 2 a. This Hamiltonian can be derived directly from Wilson’s lattice 
Lagrangian using either the transfer matrix or canonical transformation. The 
gauge fields have to satisfy the commutation relations 

[E a (x,k),U(yJ)\ = T a U(x,k)6 x J ktj , 
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[E a {x,k),U\y,j)\ = -U\x,k)T a 8 x , y 8 k , h (3) 

with T a being the fundamental representation of the ath generator of the 
gauge group. 

The wave function of the ground state is assumed to be of the form 


\Sl)=exp[R(U)]\0), (4) 

where the bare vacuum |0) is defined to be fluxless and RiU ) consists of 
gauge invariant operators such as the Wilson loops. The vacuum state with 
energy cq has to satisfy the lattice Schrodinger equation 

H\n)=e a \n), (5) 

which results in an eigenvalue equation for H 

9 9a 

Y,m, [E,. ii(f ’)]] + [£<• -ft(l r )][£,, -R(C)]} - - E Tr(V r + I'D = —ea. (6) 

l 9 p 9 

This equation can be solved by a truncation method [6, 7], in which RiU ) is 
expanded in order of graphs, 

RW) = '£R i (U), ( 7 ) 

i 

and the order is defined as the number of plaquettes involved. Denote Ri and 
Rj as the graphs of order i and j respectively, and all new graphs created by 
J2i[Eh Ri(U)] [Ei, Rj(U)] are defined as graphs of order i + j. Then the nth 
order truncated eigenvalue equation is 

n 2 

E{[Ci,[C,.E«.(C)]]+ E [E h R,m[E h R,m\ - -T.T'iu r + ul) 

l i 3 p 



( 8 ) 


The long wavelength limit of a graph is obtained by small a expansion 
of the graph. In this limit, the vacuum wave function (4) is reduced to 
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the continuum one (1). The ground state of the form (4) implies that at 
large scales, the pure gauge lattice vacuum is governed by multi monopole 
configurations. 

In SU(2) gauge group, TrU p = TrU^ and all loops with crossing can be 
transformed into loops without crossing. Then according to this rule, there 
are one graph of first order, three graphs of second order, and nine graphs of 
third order. In SU(3), however, these are no longer the case. 

In the analytical calculation, the unitary and unimodular conditions lead 
to constrains on the graphs. Any group element A of SU(3) has to satisfy 
the following condition 

^-uii Ai 2 j 2 Ai 3 j 3 tj 1 j 2 j 3 ti 1 i 2 i 37 (9) 

where a summation over the repeated indices is implied. We rewrite this 
condition as 

2 (A% = 2 (A 2 ) tJ - 2 A tJ TrA + [(TrA ) 2 - Tr(A 2 )]%, (10) 

or 

2 % = 2(A 3 ) tl - 2 (A^ijTrA + [(TrA) 2 - Tr(A 2 )]A tl7 (11) 

from which the relations among different graphs can be established. For 
example, 

2 G\ = G' 2) i ~ Cr2,l, 


6 — Gg x — 3Gg 2 + 2(73,1 


G*2,4 — ^2,3 ~ ^3,6 + G' 3 3 


Cr2,5 — Cr2,6 — Cr3,3 + G' 3 


* 2,6 


3,3 


' 3,45 


( 12 ) 

where G' and G are the graphs which can be found in Fig. 1 and Fig. 2 with 
G\ = □ = TrUp and so on. One sees that not only graphs of the same order, 
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but also graphs of different orders mix, so that the classification becomes 
rather involved. In this paper, we choose as much as possible the connected 
graphs (e.g. 6 * 21 ) as independent graphs, and try our best to transform the 
disconnected ones (e.g. G' 2 i) into the connected ones, because we think that 
the connected ones give more relevant physical information at a given order. 
The complete set of graphs up to the third order are 

R X {U) = CiGi + h.c., 

R 2 {U) = ^ 6 * 2 , 8 ' 6 * 2 ,; + h.c., 

8 = 1 

29 

Rs(U) = ^ 63 ,; 6 * 3 ,; + h.c., (13) 

8 = 1 

which graphs are plotted in Figs. 2.1, 2.2 and 2.3 respectively. Substituting 
them into ( 8 ), we obtain the 36 nonlinear equations for the coefficients 6 *i, 
C 2 ,t and C 3ti . 

Physical quantities like /i 0 and /i 2 are related to the coefficients of the 
graphs in the long wavelength limit: 


/lo — [Cl + 6(62,2 + 6*2,4 + 6*2,5) + 4(6*2,! + 6*2,3) + 15 ( 6 * 3,2 + 6*3,4 + 6*3,5) 


+ 27(6*3,7 + 6*3,12 + 6*3,13 + 6*3,14 + 6*3,19 + 6*3,20 + 6*3,29) 


+ 3(6*3,15 + 6*3,16 + 6*3,24 + 6*3,25 + 6*3,28) + 6*3,6 + 6*3,17 + 6*3,18 + 6*3,23 + 6*3, 


26 


+ 15 ( 6 * 3,8 + 6*3,10 + 6*3,11 + 6*3,21 + 6*3,22) + 9(6*3,9 + 6*3,3 + 6*3,1 + 61*3,27)] <7 4 , 


fJ - 2 — [6*3,6 — 6*3,3 + 2(63,15 + 6 * 3,16 + 6 * 3,24 + 6 * 3,25 + 6*3,28) 
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(14) 


3 

2 


(^ 3,8 + C: 


3,10 


+ ^ 3,11 


+ ^ 3,21 + (^ 3 , 22 ) 


— 36*3^ — Cs ,17 ± 26*3^8 + 63,23 — 26*3^7 + —(^ 


2,6 — ^ 2 , 3)19 


Ho and // 2 should be constants in the weak coupling limit g —> 0 or 
f3 = Q/g 2 —> 00 as required by the renormalizability of the theory. Fi¬ 
gure 3 gives a comparison between the third order results from the strong 
coupling expansion and the truncated eigenvalue equation. (The solid line is 
made by joining 100 data in the interval f3 £ [0.12,12], while the crosses are 
only the representative points.) They are consistent in the strong coupling 
region (/3 < 6 ), which implies that the calculation using the truncated eigen¬ 
value equation method is supported by the results from the strong coupling 
expansion. For larger /3, it is not surprising that the strong coupling expan¬ 
sion method no longer works. It is usually hoped that beyond the strong 
coupling region, there is a scaling region for extracting continuum informa¬ 
tion when the physical quantities become approximately constants. From 
the intermediate coupling (/3 ss 6.84) till the weak coupling (/3 ss 11.52), the 
data from the the truncated eigenvalue equation method show nice scaling 
behavior, thus suggesting the correct long wavelength continuum limit ( 1 ) of 
the vacuum wave function (4). From the results we estimate 


Ho ~ 0.5411 ± 0.0038, 


H 2 ~ -0.0781 ± 0.0024, (15) 

where the mean values are the averaged ones over the 40 data in the scaling 
region f3 £ [6.84,11.52], while the error analysis is based on the jackknife 
method (only a rough evaluation for the errors). As far as we know, the data 
for Ho and h 2 from Monte Carlo simulations or other analytic methods are 
still lacking. 

The presence of the expected scaling behavior in the intermediate and 
weak coupling regions tells us that even at finite bare coupling ^ 7 , it is possible 
to extract the continuum information along the line of constant physics, since 
it will eventually flow into the critical point g cr = 0 , i.e., the continuum limit 
when the cut-off becomes larger and larger. These results again indicate 
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that the correlation length of the continuum gauge held strength tensor in the 
ground state of QCD 3 has the order of 1/e 2 , i.e., the order of the confinement 
scale. Of course, as the continuum limit a —> 0 or equivalently f3 —> oo is 
approached, the correlation length in the lattice unit will be divergent so that 
the inclusion of higher and higher orders of the graphs is required to better 
represent the vacuum state. 

There remains the problem about the choice of set of independent gra¬ 
phs. Due to the fact that elements of the gauge group have to satisfy the 
unitary and unimodular conditions, graphs of different orders mix. There is 
ambiguity in choosing a set of independent graphs at a given order. Different 
choice of independent graphs might give different results when truncating at a 
finite order. In this paper, we have chosen as much as possible the connected 
graphs, because we think that these graphs may stand for more coherence, 
and may lead to more rapid convergence to the continuum limit. Mixing 
with lower order graphs also appeared in [5] as a result of the shifting proce¬ 
dure. Here we invoke the unimodular condition for the necessity of mixing. 
We think that mixing is essential and there might exist some certain scheme 
for most efficient approach to scaling. This problem is currently under close 
investigation. 

In summary, we have successfully applied the truncated eigenvalue equa¬ 
tion method [6] to the realistic group SU(3) with some new prescription 
scheme for the classification of the graphs, and determined nonperturbatively 
the ground state wave function. The observation of scaling for // 0 and /r 2 in¬ 
dicates the correct continuum behavior of the lattice vacuum wave function 
at relatively large scale or long wavelength. Once a reasonable form of the 
ground state is established, other physical quantities such as the spectrum 
of the excited states can be evaluated. Extension of our QCD 3 model and 
techniques to 3+1 dimensions is hopeful to yield phenomenologically relevant 
results. Such work is in progress [15]. 
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Figure Captions 


Fig. 1. Graphs G' in (12). 

Fig. 2.1. First order graph in RiU). 

Fig. 2.2. Second order graphs in RiU). 

Fig. 2.3(a). Third order graphs in RiU). 

Fig. 2.3(b). Third order graphs (continuation of Fig. 2.3(a)) in RiU). 

Fig. 3. Ho and /r 2 as a function of /3, where the triangles stand for the results 
from the strong coupling expansion, while the crosses are those from the 
third order calculation of the truncated eigenvalue equation method. 
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